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u ■ Abstract 



We present the third-order contributions to the quark-gluon and gluon-quark timelike splitting 
functions for the evolution of fragmentation functions in perturbative QCD. These quantities have 
been derived by studying physical evolution kernels for photon- and Higgs-exchange structure 
functions in deep-inelastic scattering and their counterparts in semi-inclusive annihilation, together 
with constraints from the momentum sum rule and the supersymmetric limit. For this purpose we 
have also calculated the second-order coefficient functions for one-hadron inclusive Higgs decay 
in the heavy-top limit. A numerically tolerable uncertainty remains for the quark-gluon splitting 
function, which does not affect the endpoint logarithms for small and large momentum fractions. 
We briefly discuss these limits and illustrate the numerical impact of the third-order corrections. 
Compact and accurate parametrizations are provided for all third-order timelike splitting functions. 



In this article we address the scale dependence (evolution) of the parton fragmentation distributions 
(functions) D^Hx.Q^), see Ref. IIJ for an introductory overview, at the next-to-next-to-leading 
order (NNLO) in massless perturbative QCD. Here x represents the fractional momentum of the 
final-state parton / transferred to the outgoing hadron /?; and Q^ is a (timelike) hard scale, for in- 
stance the squared momentum of the virtual photon or Z-boson in semi-inclusive electron-positron 
annihilation (SIA), e+e^ — )■ y, Z — t- /? +X, where X stands for all accessible hadronic states. 

The evolution of the fragmentation distributions is given by 

where the summation over b = q^, q^, g for z = 1 , . . . , n^ is understood, and rir denotes the number 
of effectively massless quark flavours. Unlike the functions D^{x,Q^), the 'timelike' splitting 
functions P^^ can be expanded in powers of the strong coupling tts. 

Pi {x,a,{Q^)) = a,PfJ\x) + alP^^J\x) + alpj^fix) + ... . (2) 

We normalize the expansion parameter as a, = as/{4n) and use, without loss of information, the 
standard MS scheme with the choice /j^ = /jj = Q^ for the renormalization and fragmentation 
(final-state mass factorization) scale. The system ([T]) of {2nr + l)x (2nr + l) coupled equations can 
be decomposed into 2nr + 1 scalar flavour non-singlet equations and the 2x2 flavour-singlet system 

Here ® abbreviates the convolution in Eq. ([U), and we have suppressed all functional dependences. 

The leading-order (LO) splitting functions P^^'^{x) O are identical to their 'spacelike' coun- 
terparts Ol for the evolution of the initial-state parton distributions (where the matrix in Eq. Q is 
transposed), a fact often referred to as the Gribov-Lipatov relation [4J. The next-to-leading order 
contributions P^^')^ were derived by several groups about thirty years ago Il5]-l9l|. Unlike the space- 
like case IfTOlfTTll . where the calculations can be performed via forward scattering amplitudes, the 
NNLO corrections py^)^ have eluded a direct calculation in terms of Feynman diagrams so far. 

The (three) non-singlet quantities P^s and the diagonal entries in Eq. ([3]) have been deter- 
mined by two of us a couple of years ago [[T2l[T3l via analytic continuations (see below) of the 
unfactorized partonic structure functions from the spacelike (deep-inelastic scattering, DIS) to the 
timelike SIA case, supplemented by complementary considerations based on Ref. IIT4ll . Only the 
second Mellin moments of Pgq (x) and iqg'^ (x) are fully known at this point |fT3l , since these 
quantities are fixed by the momentum sum rule and the diagonal entries. 

Beyond the leading order, there is no direct relation between the spacelike splitting functions, 
or their x —^ l/x analytic continuations (AC), and their timelike counterparts in the MS scheme. 
However such a relation exists for spacelike and timelike physical evolution kernels K{x,as) for 
photon-exchange DIS and SIA structure functions at NLO [[TSl . see also Ref. [fT6ll . Defining the 
expansion coefficients K^"' as in Eq. Q, it can schematically be written as 

AC[K^"^^{x)] ^ K^"^^{x) for n = 0, I . (4) 



In order to access all four splitting functions, we study the physical evolution kernels for the 
system Fi and F,^ of flavour-singlet DIS structure function. The former quantity is chosen, instead 
of ^2 in Ref. IITtI . since it directly corresponds to the transverse fragmentation function Ft in SIA. 
The NNLO coefficient functions for these observables are known from Refs. liT8l[T9l , see also 
Refs. Il20ll2n . F(|) is the structure function for DIS by the exchange of a scalar (|) coupling directly 
only to gluons via a term (|) G^'^G^y in the Lagrangian (such as the Higgs boson in the limit of a five 
massless flavours and a very heavy top quark), where G^^ denotes the gluon field strength tensor. 
The NNLO coefficient functions for the structure function F^^ have been calculated in Refs. lil7IJ22ll . 
while those for the corresponding fragmentation function FI are presented in Appendix A. 

The spacelike physical kernels have been discussed, for the system (F2, F^) , in detail in Ref. ifTTl . 
The timelike case is completely analogous up to a transposition of the matrices: 

^ F^ = K^®F^ = {(p^+C^®/'^)®(C^)-i)}®F^ (5) 



dlnQ^ 



with 



(") t-(") 



where c^j^' = Ci ^ = 5(1 —x) and c^- = Ci „ = 0. We have skipped the superscript T' where 
it is not needed for uniqueness in the present context. It should be noted that the normalization of 
Cj in Eq. ^ differs by a factor of 1/2 from that in Refs. II19II21II . Finally (3 in Eq. ^ is the standard 
beta function of QCD, (3 = —^oa^ + ... with Po = 1 1/3 Ca — 2/3 rir and Q = A'coiours = 3. 

For the off-diagonal entries the analytic continuation involves, besides x — )■ 1/jc and the mul- 
tiplication by a factor x due to the phase space of the detected parton in the SIA case [[T9l , a sign 
factor and a ratio of colour factors, leading to (with Cf = 4/3 in QCD) 

AC[K^''^{x)] = -^xK^^{\/x) , AC[K^^{x)] = - ^xK^'^ {\/x) . (7) 

The critical part of the analytic continuation is that of powers of In ( 1 —x), which is given by 

AC 

ln(l -x) — ^ ln(l -x) - \nx + Kin with K = or 1 . (8) 

For K = 1 the real part is taken in the end. It is not clear at all that beyond NLO Eq. ([8]) is applicable, 
in either form, to quantities such as physical evolution kernels instead of to (classes of) Feynman 
diagrams, see the discussions in Refs. ll5l[8l[T6ll. 

The NLO physical kernels for (Fi,/>) and {Ft,F^) fulfil Eq. dH) for both K = and K = 1 in 
Eq. dS]). However, for the NNLO diagonal entries we find (restricting ourselves to jc < 1) 

ACo [K[f{x)] - Ki^^ix) = 24^2poC2 i±^ \nx = 12^2|3oC^<^ W Inx , (9) 

J- -A. 

(0) (0) (2) (2)T 

and a completely analogous relation with Pgg = KL on the right-hand- side for KL and KL . 

Using K = 1 instead leads, besides ^2 Po terms with /qq (x) and Pqq (x) Inx, to the same spurious 
C,2Cp. contribution as found in the analytic continuation in Ref. [[T2|. 



There is no obvious reason why an NNLO imperfection of the AC relation should lead to an 
offset proportional to |3o and the lowest-order kernel, and why exactly the same relation should 
hold for two different kernels. Hence it seems very likely that a non-vanishing r.h.s. of Eq. ^ 
is genuine and analogous to the 'Crewther discrepancy' between the Gross-Llewellyn Smith sum 
rule in DIS and the Adler function in e+e annihilation observed and discussed in Ref. Il23l . see 
also the all-order proof and the recent explicit fourth-order calculation in Refs. Il24l . 

Therefore, if one tried to fix the so far unknown ^2 terms of the off-diagonal splitting func- 

(2) (2) 

tions by imposing Eq. (H)) for K^J , Kl^ and their timelike counterparts, one should find an offset 
proportional to Po in the known second moments which is then compensated by right-hand- sides 
analogous to that of Eq. ^. Carrying out these calculations indeed leads to 

'ACo[Kif]-4^^]^^^ = ^2|3oC^(-^Q-16C^) , (10) 

ACo[K^f]-Kll^]^^^ = C2|3o«/(9Q-fC^+4|3o) . (11) 

For the required generalization of these relations to all values of N, we consider also the K = 1 
continuation which, while again leading to spurious non-Po terms, appears to provide the right 
correction terms for Eq. (flOl) . 

ACo [K^l) -4^^ = -6C2|3o/'iJ^(x) [2Q(1 -Inx) + C^] , (12) 

as well as the |3q part of Eq. (fTTI) . We assume that also the first two terms on the r.h.s. of Eq. (fTTI) 

correspond a to combination of iqg and Pqg \nx. (Poly-)logarithms with higher weight cannot 
occur since this is an n. contribution and thus restricted to an overall weight of three. Imposing 
also other constraints discussed below we arrive at 

ACo[Klf]-Kll^ = C2|3o<^(^)[8(Q-C^)-121nx(Q-2Q) +6po] • (13) 

Also the remaining uncertainty of the coefficients of \nx in this relation will be addressed below. 

We are now ready to present the (only marginally provisional) results for the NNLO timelike 
splitting functions. For completeness we first recall the corresponding LO and NLO results: 

p4°^^W = ^qg^^W = 2n^pqg(x) = 2n^ (1-2^ + 2x2), (14) 

Pij)^(z) = /'iJ)^(x) = 2CpPgq{x) = lCp{2x-'-2 + x) (15) 

and 

P^^'^{x) = ACpHf ( - 6 + 23/2x - 10x2 - (5/2 -2x- 2x^) Hq - (1 - 2jc + Ax^) Ho,o + 2Hi 
+ 2pqg(x)(-C2 + Hi/2-3Hi,o-Hij+H2))-4/3n/(2 + 2pqg(x)(2/3 + Ho-Hi)) 

+ 4Qn^(-(20jc"^-13 + 95x-178jc2)/9-4^2^-(4 + 34x-f4jc2)/3Ho 

+ (2-M2.^)Ho,o-2Hi-277qg(-.^)H_i,o + 2;.qg(x)(2Hi,o-5/6Hi-FHi,i-2H2)), (16) 

P^^'^ix) = 4C^((9x-l)/2-(8-x/2)Ho + (2-x)Ho,o-2xHi+2pgq(x)(2Hi,o + Hij 
- 2H2)) + 4C^Q (+ 17/9x^1 + 5 -X- 44/9x2 +4^2 -(6x"^- 8 -9x- 8/3^2) Ho 
-(8x-i+4 + 6x)Ho,o + 2xHi-2;?gq(-x)H_i,o-2;?gq(x)(3Hi,o + Hi,i-H2)) . (17) 



//,7, are the harmonic poly logarithms (HPLs) as defined in Ref. [f25l| with i/o,i,o,i(-'f) — ^2,2 etc. 
Our new third-order contributions to Eq. (O are given by 

pS'iA = 

CJHfi - 227/4 - 1967/2X + 1069;c2 - (288 - 352;c + 160;c^)Hi,o,o + (34 - 32x - Ax^)ni^i 

-(240-64jc + 256jc2)H_2,o-(180-192x+144x2)Hi,i,o-(176-288x + 416x2)H2,o,o 

-(124-64jc + 48x2)Hi,2-(120-240;c + 288;c2)H2,i,o-(104-208x + 224;c2)H2,2 

-(96-64x + 256x2)H_3,o-(76-64jc + 48x2)Hij,i-(56-112x+128jc2)H2,i,i 

_ (40- 144JC+ 160jc2)Ho,oC2- (28 + 144jc- lex^)^^ - (2- 116jc + 324x2)H2 

+ (156 -816x- 640x^)^3 + (12 + 272x-176x2)HoC2 + (16 -256x+128x2)Ho,o,o,o 

+ (32 + 128jc + 256jc2)H_2,o,o + (36 - 584x - 192x2)Ho,o,o + (64 + 320x + 256x2)H_i,o,o 

+ (40 -304x+160x2)HoC3 + (56 -1904/5X + 352/5x^)^1 + (64 -128x+192x2)H3,i 

+ (64 + 64x+128;c2)(H_2C2 + 2H_2,-i,o)-(496 + 240x-224x2)H_i,o 

+ (92- 1143X+ 110x2)Ho + (96- 192x + 352x2)H3,o + (104- 144jc + 224jc2)H2C2 

+ (104 - 80jc+ 16x2)H2j + (106 - 564x + 564;c2)^2 + (107 - 42x- 76x2)Ho,o 

+ (108-64x-80x2)HiC2 + (128 + 160x + 32x2)(H_i^2 + 2H_i_i,o) 

+ (180 - 480x + 318x2)Hi + (184 - 128jc - 1 12x2)H2,o + (294 - 560x + 468x2)Hi,o 

+ (8-16x + 32x2)H4-96pqg(-x)(H_i^_2,o + 2H_i^_i^o,o + H_i,oC2-ll/6H_i^OAo) 
+ 1677qg(x) (3IH1C3 - 6Hi,_2,o - 7Hi,oC2 + 21Hi,oao - Hi,iC2 + Hi,i,o,o + 3Hi j,i,o 
+ 2Hi,i,i,i+5Hi,i,2 + 19Hi,2,o + 4Hi,2,i+5Hi,3)) 
+ C4C^n^(17597/36 - 220/3;c-i + 2659/6jc - 3092/3^^ - (288 + 448jc + 64jc2)H_i _i,o 

- (592/3 - 1352/3JC+ 192jc2)H2,i - (2422/27 -30203/27jc+ 17918/27x2)Ho 

- (188-32x^1 -272x+ 120x2)HiC2 - (1007/9 + 6254/9x-2444/9x2)Ho,o 
-(160 + 160x + 320x2)H3,o-(176-192x + 384x2)H3,i-(96 + 160x-32;c2)H_iC2 

- (64 + 32x+ 128^2) (H_2C2 + 2H_2 -1,0) + (28/3 + 2032/3x + 280/3x2)H2,o 

- (488/9 - 304/9x^1 - 868/9jc+ 1904/9;c2)Hi ^ - (148/3 - 1832/3x- 640/3x2)Ho,o,o 

- (48 + 64x + 64x2)H_i,2 - (32 + 64x)Ho,o,o,o - (16 + 96x"* + 504jc - 368x2)Hi,o,o 

- (44/5 -2632/5X+ 128x^)^1- (56 -16x+160x2)HoC3 + ( 16 -80x-144x2)H_i,o,o 
+ (32 -256jc + 256jc2)Ho,oC2 + (64 + 96x+128x2)H4 + (28/9 + 4540/9JC- 2620/3x^)^2 
+ (268/3 - 32jc"1 - 608/3x + 728/3x2)Hi,2 + (-260/3 - 176/3x - 8x2)HoC2 

+ (2788/27 - 80/27x^1 + 19660/27x - 19970/27x2)Hi + (388/3 + 256/3x - 40/3x2)H3 

+ (136-1040x + 672x2)H2,o,o + (136-368x + 320x2)H2,ij-(200-592x + 512x2)H2C2 
+ (424/3 -32/3x-i-536/3x + 544/3x2)Hi,i,i + (2204/9 -3772/9x + 3176/3x2)H2 
+ (452/3- 32x^1 -1024/3x+1048/3x2)Hi,i,o + (168 -624x + 448x2)H2,2 
+ (184-656x + 512x2)H2,i,o + (144 + 32x + 384x2)H_3,o + (336 + 16x+160x2)H_2,o 
+ (1340/3 + 64x^1 +4904/3x + 488x2)^3- (212- 304/3x^1 -2036/3x + 3284/3x2)Hi,o 
+ (672 + 544x - 16x2)H_i,o - (96x + 128x2)H_2,o,o - 8pqg(-x) (10H_i,_2,o - 8H_i,_i,o,o 



- 10H_i _i^2 + 5H_i^3 — 12H_i _i _i^o + 4H_i _i^2 + 3H_i^o,o,o — 12H_i^2,o — 4H_i_2,i) 

- 8pqg W (99Hi^3 - 10Hi,_2,o - 10Hi,oC2 + 7Hi,o,o,o - 4HijC2 - 26Hi,i,o,o + 2Hi,i,i,o 
+ 8Hiaaa+6Hi,i,2 + 30Hi,2,o + 4Hi,2,i-14Hi3)) 

+ ci«^(2057/9 + 2092/81jc"^ - 672JC + 53753/81x2- (196 + 2416/3x-328/3jc2)^3 

- (3940/27 - 152/21x-^ + 8416/27x - 7448/27;c2)Hi + (288x - 96^2) (H2,i,o + H2,2) 

- (1262/9 -112/9jc"^-4828/9jc + 2900/3jc2)H2- (100 -64/3jc~^ + 312x-72x2)H2,o 
-(352/3 + 64/3x^1 +704/3x + 80x2)H_i,2-(100 + 776x-96x2)H4- (64+ 128;c2)H_2,2 

- (2996/27 - 40jc"1 + 4832/27x - 5044/9jc2)Ho - ( 104 - 432jc + 32x2)H_2,o,o 

- (80-256JC+ 192/)H2,i,i - (598/9+ 128/9x"' +3508/9x + 7208/9x2)Ho,o 

- (196/3 -32/3x~^-344/3x + 400/3x2)Hi,i,i + (168 -432x + 384x2)H2C2 
-(310/3- 544/9x^1 -308jc + 536/3x2)^2- (56 -496jc)H_3,o + (64 -1264x)Ho,o,o,o 
-( 152/3 + 32/3jc"^+256/3x + 248x2)H_i,o- (28 -632x + 96x2)Ho,oC2 

+ (96 + 64/3x"^-484/3x + 824/3x2)HoC2 + (134/9-128/9x"^-700/9x + 652/3x2)Hij 

+ ( 16 - 352jc - 64x2)H_2,-i,o + (V28/3 + 224/3jc"1 + 188/3x - 304/3x2)Hi,o,o 

+ (64/3 - 128/3x^1 +2680/3X- 192x2)Ho,o,o + (64/3 - 64/3x"' +292x+ 120x2)H3 

+ (34 + 252/5x + 704/5x2)^2^ (160/3 + 32x"^+304/3x-472/3x2)(Hi, 1,0 + Hi,2) 

+ (184/3 + 64/3x"^-344/3x + 832/3x2)H_2,o + (64 -96x+160x2)H3,o 

+ (440/3 - 64/3x^1 - 664/3x + 464/3x2)Hi ^2 + (72 - 176x + 96x2)H_2C2 

+ (92 + 408x - 32x2)H2,o,o + (320/3 - 64/3x"i + 1 12/3x - 208x2)H_i,_i,o 

+ (128 + 32/3x"^-400x + 640/3x2)H2,i + (16 + 960x)HoC3 + (176-192x + 320x2)H3,i 

+ (512/3 + 32/3x^1 +760/3x-24x2)H_i^2- (562/9 + 8/9x^1 -92/9x-704/3x2)Hi,o 

-(872/3- 128/3x^1 + 520/3x-304/3x2)H_i,o,o + 8/?qg(-x)(17H_i^3 + 14H_i_2,o 

-26H_i,_i^2-12H_i,_i,_i,o + 8H_i,_i,o,o + 20H_i,_i,2 + 14H_LoC2+H_i,oAO 

- 8H_i,2,o - 4H_i,2,i) + 8pqg(x)(31HiC3 - 2Hi,_2,o + 6Hi,oC2 - 15Hi,o,o,o - 2Hi,iC2 
-20Hia,o,o-4Hi,i,i,o + 4Hi,i,i,i -4Hia,2-4Hi,2,i - 16Hi,3)) 

+ C^«/( - 4847/54 + 200/27x"^ -2375/27X + 4066/27x2 - (32 - 64x)Ho,oAO 

- (416/9 -376/9x + 88/3x2)H2 + (1684/27 + 512/27x"^-4022/27x-4480/27x2) Ho 

- (808/27 - 560/27x + 248/27x2)Hi + (902/9 + 128/9x^1 +956/9x+ 1400/9x2)Ho,o 

- (64/3 - 128/3X+ 176/3x2)H3 + (272/9 - 64/9x"i + 8/9x + 32x^)^2 

+ (48- 64/9x^1 -496/9x2)H_i,o + (344/9 -376/9x + 328/9x2)Hia 

+ (80-160x+176x2)/3H2,o-(32-112x-96x2)/3Ho,o,o + (64-128x+144x2)/3H2a 
+ (32 -64x + 64/3x2)H_2,o + (80 -352x+144x2)/3HoC2 + (304/3 -608/3X+ 160x^)^3 

+ (192-248x + 232x2)/3Hi,o-8/3pqg(x)(3Hi^2 + 5Hi,2 + 7Hia,o + 5Hiaa-12Hi,o,o)) 

+ C4n/(- 14/9 -44/9x"'-1216/9x + 916/9x2- (40 -352/3X + 272/3x^)^3 

- (296/9 + 64/9x^1 + 416/9X + 536/9x2)Ho,o - (32 - 96x + 224/3x2)H2a 

- (64/3 + 640/3x)Ho,o,o + (44/27 - 256/27x^1 + 752/27x - 2212/27x2)Ho 



+ (100/3 + 32/9jc"1+208/3x-16/9jc2)^2- (24- 112/3jc+112/3jc2)HoC2 

- (320/9 - 32/9jc"1 -544/9x + 200/3x2)Hi,i - (64/3 - 128/3jc + 64/3jc2)H_2,o 
+ ( 16/9 + 32/3x"^+496/9jc-664/9jc2)Hi,o + (356/9 -1360/9x + 632/9jc2)H2 
+ (592 -320^-^ - 1232x + 2632jc2)/27Hi + (104/3 - 16jc+ 176/3jc2)H3 

-(144-32jc"i-96jc-512jc2)/9H_i,o + (8 + 80jc-16jc2)H2,o + 32/3pqg(-Jc)(2H_i^2 
+ 2H_i,_i,o + H_i,o,o-H-i,2)-8/3pqg(x)(4Hi^2 + 19Hi,o,o + 2Hi,i,o-5Hi,i,i+2Hi,2)) 

+ «/(16/9 + 8/3[Ho-Hi+;?qg(x)(C2-l+2/3(Ho-Hi)+Ho,o-Hi,o + Hij-H2)]) (18) 
and 

Cj{ - 1915/2 + 794x^1 +731/4X- (464- 320jc"^ -200;c)H2C2 - (290 + 251x)Ho,o 

- (352 - 144jc"1 - 252x)Hi^2 - (264 - 180x"i - 140x)Hi,o - (240 - 256jc"1 - 120x)H3,i 

- (230-228jc"1 - 16jc)Hi - (224 + 256x^1 - 148jc)^2 - (176- 128x"1 - 88x)H3,o 

- (776/5 - 1216/5x^1 - 84/5x)C| - (128 - 64x)H4 - (104 - 40x)H2,o + (248 - 28x)H3 

- (96 + 720x^1 + 196x)^3 - (80 + 256x^1 -208x)H_i,o- (80- 128x"i +200x)HoC3 

- (72 - 60x^1 - 80x)Hij - (72 - 56x)H2,i - (64 + 32x)(H_iC2 + 2H_i _i,o) 

- (64 + 96x)H_3,o - (64 + 64x)Ho,oC2 - (64 + 128x"i + 64x)H_2,o,o + (208 - 140x)H2 
-(32-48x-i-4x)Hi,i,o-(32-96x)Ho,o,o,o-(24+192x-i+20x)HoC2 

+ (197 + 208x^1 +29x)Ho + (16 -52x)Hijj + (128 + 128x^1 +96x)(H_2C2 + 2H_2-i,o) 
+ (128 - 192x^1 - 1 12x)H_2,o + (128 - 48x"i - 60x)Hi,2 + (128 + 48x"i - 36x)Hi,o,o 
+ (80 - 64x^1 - 40x)H2,i,o + (144 - 64x'i - 72x)H2,2 + (160 + 40x)Ho,o,o 

- (16 - 128x^1 -40x)H2,o,o + (208 - 192x"i - 104x)H2,i,i + (224 + 96x^1 +96x)H_i,o,o 
+ 16pgq(-x)(2H_i,_2,o + 4H_i,_i,o,o-H_i,oAo) + 16pgq(x)(2Hi^3-2Hi,_2,o + 4Hi,oC2 
-6Hi,o,o,o + 3Hi,i^2-6Hi,i,o,o + Hij,i,o + 2Hi, 1,1,1 -3Hi, 1,2 + Hi,2,o-4Hi,2,i-8Hi,3)) 

+ C4C^(1735/6 + 200x^1 -481 l/12x- 140/3x2 -(336 -448x^1 -312x)H2,2 

- (2338/3 - 80x^1 +441x+ 1408/9x2)Ho,o - (1760/3 - 336x"i - 832/3x + 64/3x2)Hi,o,o 

- (336 - 384x^1 - 216x)H2,i,o - (8980/9 + 5800/9x^1 +4576/9x + 352/3x2)H2 

- (976/3 - 904/3x^1 - 404/3x + 32x2)Hi,2 - (320 - 256x"i - 208x)H2,o,o 

- (272 + 832x^1 + 376x)HoC3 - (272 - 320x"i - 184x)H2,i,i - (208 + 448x"i + 344x)H3,o 

- (192 + 48x^1 +72x)H_i,o,o - (144-96x^1 - 128x)H_2,o - (128 - 32x"i + 16x)H_2C2 

- (128 + 64x^1 +96x)H_2 -1,0 -(256/3- 88x^1 -332/3x + 32/3x2)Hi, 1,0 

- (608/5 + 3664/5x^1 + 2144/5x)C| - (88 - 192x"i + 268x)H_i,o - (32 - 48x)H_3,o 

- (152/3 - 160/3x^1 -232/3x + 32/3x2)Hi, 1,1 + (552 + 264x^1 +424x+ 128/3x^)^3 

+ (136/3 + 640/3x^1 - 680/3x - 32x2)H3 - (160/9 - 848/9x"i + 280/9x + 352/9x2)Hi,i 

- (64 + 24x^1 +32x)H_i^2- (16 + 128x^1 +40x)Ho,oC2- (16 + 320x^1 +216x)H3,i 
+ (72 + 192x^1 + 128x - 32/3x2)HoC2 + (368/3 - 268/3x + 64/3x2)Ho,o,o 

+ (424/3 -464/3x^1 +28/3x + 32x2)H2,i + (600- 1280/3x^1 +96x + 224/3x2)H2,o 



+ (488/3 + 620/9jc"1 - 200x - 704/9jc2)Hi,o + (192 + 48jc"1 + 128x)H„i,_i,o 

+ (272 -608jc"1-456x)H2C2 + (336 + 512x^1 +376jc)H4 + (160 + 48x"^+96jc)H_i,2 

- (32 + 32jc) (H_2,o,o - 2Ho,o,o,o) + (40301/27 + 244x"' + 1891 l/27;c + 8932/27x2)Ho 

+ (672- 1760/3X-1 - 180jc + 224/3x2)Hi^2 + (3452/3 + 528x"' +436/3;c+ 1760/9jc2)^2 
+ (26008/27 - 39092/27jc~i + 3988/27jc + 5908/27jc2)Hi + 8pgq(-x)(4H_i,3 + 1 lH_i^3 

- 18H_i _i^2 — 8H_i^2,o — 4H_i^2,i — 4H_i _i^o,o+ 12H_i _i^2 + 12H_i^oC2 + H_i^o,o,o 
+ 10H_i,_2,o - 12H_i,_i,_i,o) - 877gq(x)(4H_2,2 + SlHi^g - 10Hi,i,o,o - 14Hi,3 - 4Hi,2,i 
+ 3Hi,oAO + 12Hi,i^2 - 6Hi,_2,o - 2HkoC2 + 8Hi,i,i,i + 2Hi j,2 + 34Hi,2,o + 14Hi,i j^o)) 

+ CiQ(- (12 + 752/3x^1 -532/3X+ 128/3x2)HoC2 - (152+ 128x^1 + 92x)H4 

- (22270/27 -31504/27x^1 +688/27X+ 1924/9x2)Hi - (144 + 64x^1 +8x)H_2C2 
-(17798/27 -3514/9x^1 +3436/27x-4136/9x2)- (400- 160x-1-148x + 640/3jc2)^3 
-(576- 1504/3x^1 +424x + 512/3x2)H2,o- (1292/3 + 32/3x^1 +284/3x+128/3x2)H3 

- (848/3 - 1856/3x^1 + 3344/3x + 224/3x2)Ho,o,o + (572/5 + 1984/5x^1 + 314x)C| 

- (32612/27 - 13384/27x^1 +23630/27X+ 12964/27x2)Ho- (160-48x)H_2 -i,o 
-(168 + 176/3x^1 + 160x + 32/3x2)H_i^2- (280-344x^1 +32x + 64x2)Hi^2 

- (848 + 5176/9x^1 +2800/9x + 352/9x2)^2 + (152- 128x^1 -156x)Ho,oC2 

-( 1280/9 + 248/3x^1 -1052/9x-704/9x2)H_Ko- (128 -5 12x^1 -632x)Ho,o,o,o 
-(328/3 + 272x^1 -544/3x+128/3x2)H_i,o,o + (128- 192x^1 - 112x)H2,i,o 
+ (384 + 576x^1 +480x)(H3,o + l/3H3j) - (296/3 -448/3x"i +260/3x+ 128/3x2)H2,i 
-(280/3 + 752/3x^1 -928/3x + 256/3x2)H_2,o + (64- 128x^1 -80x)H2,i,i 

- (16 - 352/3x^1 + 160x - 64/3x2) H_i_i,o + (1300/3 - 1024/3x^1 - 1004/3x)Hi,o,o 
+ ( 104/3- 160/3x^1 -76/3x + 32/3x2)Hi,i,i + (48 + 320x^1 + 248x)H2C2 

+ (11848/9 + 3920/9x^1 +7280/9x + 3688/9x2)Ho,o + (112- 256x^1 -248x)H_3,o 

+ (368/3 - 424/3x^1 - 352/3x + 32/3x2)Hi j,o + (160 + 352/3x"i + 80x + 64/3x2)H_i,2 

+ (240 - 96x^1 - 136x)H_2,o,o + (1 160/9 - 4136/9x^1 - 604/9x + 760/3x2)Hi,o 

+ (64 - 256x^1 - 176x)H2,2 + (352 + 416x^1 + 240x)HoC3 + (232 - 160x"i - 76x)H2,o,o 

+ ( 192- 248x^1 -72x + 32x2)Hi,2 + (6896/9 + 4664/9x^1 +4220/9x+1552/9x2)H2 

+ (808/9 - 532/3x^1 - 440/9x + 560/9x2)Hi,i - Spgq{-x) (4H_2,2 - H_i^3 + 22H_i _2,o 

- 2H_i _i^2 — 12H_i _i _i^o + 12H_i _i^o,o — 4H_i _i^2 + 10H_i^oC2 — 21H_i^o,o,o 

- 12H_i,2,o - 4H_i,2,i + 4H_i,3) + 8pgq(x)(41Hi^3 - 6Hi,_2,o - 8Hi,oC2 + 35Hi,o,o,o 
+ 6Hi,iC2+10Hi, 1,0,0 + 12Hi,i,i,o + 4Hi,i,i,i+8Hi,i,2 + 40Hi,2,o + 4Hi,2,i + 10Hi,3)) 

+ C|r^( -12803/27 + 14008/81x^1 +30475/54X- 21784/81x2 + (64 -32x)Ho,oAO 

- (204+128/3x^1 - 142/3x + 256/9x2)Ho,o- (176/3 + 128/3x"i - 160/3x)Ho,o,o 

- (484/27 - 448/27x^1 - 464/27x)Hi + (160/9 - 160/9x^1 - 104/9x)(Hi,i - 2H2) 

- (2330/27 - 200/9x^1 - 2494/27x - 3040/27x2)Ho - 8/3pgq(x) (3Hi^2 + 5Hi,o 
+ 5Hi,o,o + Hi,i,o-Hi,i,i+Hi,2 + 2H2,i-4H3)) 



+ QC^n^ (332/27 + 3778/8 1x"1-2624/27jc+1 892/8 1jc2 + (64 -80jc"1-40jc)C3 

- (256/9 - 656/9jc"1 +464/9X+ 128/9x2)Ho,o- (160/9 + 160/9jc"1 + 104/9x)H_i,o 

- (160/9 - 160/9x^1 - 104/9;c)Hi,i + (32/3 + 128/3x"i + 128/3x)Ho,o,o 
+ (136/9 - 160/9x^1 - 104/9x)H2 + (448/27 - 448/27x"^ - 392/27x)Hi 

+ ( 1280/27 + 1156/27jc~i+308/27x + 64/9x2)Ho- (320/9 -320/9x-1-136/9jc)Hi,o 

- (32/3^-1 + 16/3x)HoC2 + (8/3 - 160/9x-i - 176/9^)^2 + 16/3pgq(-x) (H_2,o 

+ H-i,o,o) + 8/3pgq(x) (6Hi^2 + 10Hi,o,o + 2Hi,i,o - Hi,i,i +H2,i - 2H3)) . (19) 

Unlike their spacelike counterparts, the off-diagonal timelike splitting functions show a double- 
logarithmic enhancement of higher-order terms not only for x — t- 1 , but also for x — )■ [261 with 

Pq^g'^^W = ^C^pnfLl + [f Q^ + f (Cf-n^)] n^U + 0(1), (20) 

Pg^''{x) = -4Q^C^Lf-8Q^QLi +0(1) (21) 

and 

xPi,^'' (x) = - f Qn^ + {xLD , (22) 

x/'g^q^^(x) = -16C^QL^-24C^QLo + f C^Q + 0(x4) . (23) 

Here and below we use the abbreviations Li = — Hi (x) = In ( 1 — x), Lq = Ho(x) = Inx and Q^ = 
Ca—Cf. The large-x and small-x limits of the new result (IT8]) and (fT9l ) are given by 

■Pqg {x) = jQ^nyLj + [-g-Q/7 + -9-Q/7(C^— «y)J Wy^l 

+ [(63i-8C2)ci^ + (^-16C2)QFCF-^QF'^^ + f(C^-n^)2]n^Lf 
+ [(^--¥^2 + 16C3)ci^ + (^-84C2 + 96C3)Q^C^ + 421c2 (24) 

-{W-f^2)CAFnf-'-fC,nf+fn}-fl^2CpiCf-nf)]nfL,+0{\), 

Pgq [x) = 3 CpC^f Lj + [ -g- Cap — g Cap [Cp — rij:) J C^ Lj 

+ [fCApnf-{^-f-%^2)Clp-{'-f-l6^2)CApCp]CpL\ 
+ [(-^ + fC2 + 80C3)ci^-(^--fC2 + 48C3)Q^C^ (25) 

+ (W-f ^2)QF^^+f(l+2^2)C^(C^-n^)]C^^i +0(1), (26) 

x/qg (x) = gCAfljrLQ — [ -g- CAflj- + -5- (C4 — 2C^) Hy J Lq 

+ [(40 + f C2)QV-|f (Q-2C^)n2] Lo + 0(1) , (27) 

xPgq (x) = — CpCA Lq + [-g- CpCA + -9- Cpnj:{CA — Cp)\ Lq 

+ [(40-64^2)c2Q + (^-64C2)C^Ci-fc2n^ + ^C^Qn^]L2 
+ [(^-¥^2 + 416C3)C^Ci + (244 + 192C2-832C3)C|Q (28) 

+ (208-192^2 + 128^3)C^ + (^-fC2)C^Qn/ + ^c2n^] ^0 + 0(1). 
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The leading logarithms in all four large-x limits above are identical to those of the correspond- 
ing spacelike splitting functions, in agreement with the all-order prediction in Ref. Il27l . As in the 
spacelike case ifTTIl . all double-logarithmic contributions to these equations vanish for C^ = Cp. 
Only one coefficient, that of In' i}—^) in Eq. (l24l) . does not vanish in the supersymmetric limit 
CA=Cf = Hf. Also this feature is analogous to the spacelike case discussed in Ref. ifTTTl . Eq. (l25l) 
is in complete agreement with the corresponding result of Ref. Il28ll . see also Ref. ||29ll . This would 
not be the case if a term with In ( 1 —x) were present on the right-hand-side of Eq. (fT2l) . 

The coefficients of the leading small-x logarithms of Pqg and P^„ at NLO and NNLO are larger 
and smaller, respectively, by a factor Ca/Cf than those of P^^ and Pgg in Eqs. (13) - (15) of 
Ref. lfT3l . For Pj and PI„ this relation and the corresponding coefficients have been derived to all 
orders in Refs. [|26ll . Unlike at NLO, P^„ and P^„ are suppressed by only one power of ln(l— jc) 
relative to Pj„ and PI„ at NNLO, and presumably some or all higher orders. 



We now turn to the additional constraints mentioned above Eq. (1131) . These are provided by the 
well-known supersymmetric relations for the choice Ca = Cf = rir of the colour factors leading 
to a fA^ = 1 supersymmetric theory. On the one hand, we can investigate the combinations 

a;;:^ W = /'cSf W ^pI^\x)-pI^^\x) -pI^\x) with a = 5, r (29) 

which vanish at LO, while at NLO and x <\ these quantities are given by 

^^^\x) = |jc-i + (4-8;c-16;c2)ln;c + f -f ;c + 28;c2, 

^^^\x) = -|;c^' -(4-8x-16x2)lnx + f + f x-4x2 (30) 

in the MS scheme lIH, see also Ref. [f30l . Obviously Ag and A^ are much simpler than the indi- 
vidual NLO splitting functions. A further simplification is obtained by adding these two quantities, 

^^l\x) + ^^■l\x) = 12-24x + 24x2 = 12pqg(jc) , (31) 

where Pqg{x) has been defined in Eq. (fT6l ). Using the results of Refs. lfT0l - [T3l and Eqs. (flSl ) and 
( fT9l ) the corresponding NNLO quantities are found to be (also at x < 1) 



Af[x) = _21n3x-91n2jc- (fx"' + ^ + 24^2)lnJC+ ... + 81n(l-jc) (32) 

^^^\x) = -2\n^x+{\6x-^-2\)\r?x+[^x-^ + ^ + 2A(,2)\rix+ ... + %\n{\-x) , 
where we have suppressed all contributions which are regular for x— t-O and x— )■ 1 for brevity, and 

A^s^^ {x) + A^p^^ {x) = - 24 ^2 Pqg (x) + non-^2 terms . (33) 



The latter means the absence of C,2 in the expansion about jc = 1 to all orders in ( 1 —x) , cf. Ref. ll3TI . 

The absence of C,2X^^lnx also in the second line of Eq. (|32|) provides a check of the coeffi- 
cient of Injc, and hence (due to the second-moment constraint (flOl) ) of the whole Ca coefficient 
in Eq. (fT2l) . An additional \nx term in Eq. (fT3] ). except with a prefactor Ca—Cf, would spoil the 
symmetry between the two lines of Eq. (|32l ) and hence also Eq. (|33] ). Obviously non-(C/7— Ca) 
terms not proportional to Pgq{x) on the r.h.s. of Eq. (fT3l) would also conflict with the latter relation. 



A second aspect concerns the analytic structure of the physical evolution kernels, in particular 
the differences between the analytic continuations (x — )■ l/x) of the spacelike and timelike ones. 
From Eq. @ we define the matrix A^^") as 

AK^"\x) = AC[if(«)^(x)] -i^W^W , (34) 

with the restriction x < I and entries according to Eqs. ©, ©, (fT2l) and (fT3l) . It is interesting 
to investigate whether Eq. (|34|) directly respects Gribov-Lipatov reciprocity. In Mellin space, this 
feature implies that the corresponding expressions are functions only of the product A^(A^ + 1) of 
the Mellin variable A^, i.e., parity preserving, a fact already exploited in the large-x (large-//) limit 
of Ref. Il28l . The eigenvalues A,,, / = 1,2 of A^^") are determined from the characteristic equation 

lj-litr:{AK)+det{AK) = 0, (35) 

and, following ||29l , it is sufficient to study the conditions which Gribov-Lipatov reciprocity im- 
poses on the trace and determinant of A^^"^ in Eq. (I34l) . i.e.. 



AC 
AC 



tr(AK^"\x)) - tr(AK^"\x)) = 0, (36) 

det(AK^"\x)) -detfAK^"\x)) = 0. (37) 



In QCD, Eq. ^ is fulfilled to NNLO due to Eq. Q, while Eq. ^ is not. However, in the 
supersymmetric limit, Ca = Cf = fir, also Eq. (ITT]) holds to NNLO as a result of non-trivial 
relations between the coefficients in Eqs. (fT2l) and (fT3]) . Eq. (ITTl) provides thus a further constraint 

(2)T (2)T 

on iqg and igq , again except for contributions proportional to Ca —Cf, which vanish trivially 
in the transition to a rA(; = 1 supersymmetric theory. 

In summary, these considerations are still not sufficient to definitely fix the right-hand- side of 
Eq. (fT3l) . As an estimate of the remaining uncertainty we suggest to use the offset 

5/'4'^^(x) = ±2i;2MCA-Cf)in+24\nx)Pi^^{x) . (38) 



The functions (1181) and (1191 ), the former including the error band due to Eq. (|381) . are shown and 
compared to their spacelike counterparts in Fig. 1; and the LO, NLO and NNLO approximations 
to Eq. ([2l) are illustrated in Fig. 2 at the typical scale Q^ ~ M|. As in the diagonal cases lfT3l . the 
higher-order corrections are much larger at small values of x than in the spacelike case. The small-x 

(2]T 12)T (2]T (2]T 

behaviour of Pqg and Pgq^ is similar to that of /qq and Pgg , respectively, with particularly 

(2)T (2)T 

large cancellations between the powers of \nx occurring in igq and igg . 

For the use of the NNLO splitting functions in numerical analyses we have prepared, analogous 
to Refs. ifTOinn . compact and accurate parametrizations not only of the present results, but also 
of the non-singlet and diagonal quantities derived in Refs. lfT2l[T3l . These parametrizations can 
be found in Appendix B. Corresponding Fortran files, and Form files of our main results, can 
be obtained by downloading the source of this paper from the arXiv servers or from the authors 
upon request. This includes the rather lengthy (even or odd) integer-A/^ Mellin-space expressions 
in terms of harmonic sums ll33l which we have not presented here for brevity. 
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Figure 1 : The third-order timelike off-diagonal (quark-gluon and gluon-quark) splitting functions 
for five flavours, multiplied by x and divided by 2000 ~ (47r)^ for display purposes. The remaining 
uncertainty of the former quantity is indicated by the dash-dotted curves. Also shown are the 
respective leading small-jc contributions and the corresponding spacelike splitting functions. 
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Figure 2: The resulting perturbative expansion of the timelike quark-gluon and gluon-quark split- 
ting functions, again multiplied by x, at a typical value of the strong coupling constant. 
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To summarize, we have performed an indirect determination of the hitherto unknown off- 
diagonal NNLO timelike splitting functions (fTSi) and (fT9l) for the evolution of parton-to-hadron 
fragmentation functions D^Hx, Q^). We expect the remaining uncertainty of the former quantity, 
estimated in Eq. (l38l) and illustrated in Figs. 1 and 2, to be phenomenologically acceptable. Hence, 
combining these results with those of Refs. lfT2llT3llT9l , NNLO analyses are now possible of data 
on the transverse fragmentation function (but not yet its longitudinal counterpart, where also the 
third-order coefficient functions are required) in semi-inclusive electron-positron annihilation. 

{2]T 

The remaining uncertainty of /qg does not affect the logarithmically enhanced large-jc and 
small-x contributions. We expect that these results can be useful to improve the corresponding 
resummations. In fact, an extension of the generalized large-x resummation of Ref. ||32]| to the 
present timelike case will be presented in a forthcoming publication. 

{2)T (2)T 

Further research is required to completely fix /qg and to check our result for igq . A direct 
calculation of the leading-/? ^ contribution to the former quantity should be possible, but would not 
address the critical contributions, while a full jc-dependent diagram calculation beyond these terms 
appears to remain formidable task. A computation of the N = 4 and N = 6 moments of the nj and 
rir parts, respectively, of Eqs. (fTSi) and (fT9l) presumably would be sufficient and may be feasible, 
e.g., generalizing the approach used in Ref. [[2T]|. in the foreseeable future. 
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Appendix A 

The first- and second-order coefficient functions for SIA via an intermediate scalar in Eq. ^ read 

cl'^ix) = Q(ll/3(l+x + x2-xfi)+pggW(8Ho-4Hi) + (67/9 + 8C2)5(l-x)) 

- nf{2/3 (1 +x + x^-x{^) + 10/95(1 -x)) , (A.l) 

4!qW = «/(4^-V^'+;^qgW(4Ho-2Hi)) (A.2) 

withpgg(x) =x^^ —2+x — x^ +x^^^, where /(x)xf ^ = f{x)/{l—x) has to be read as a -I— distri- 
bution if /(x) does not vanish at x = 1, and 

c^^lix) = Cj ((5099x^1 -3301jc2-2570jcfi)/27- 46/3 + 1858/9JC- (140/3 + 32/3X 
+ 44/3 (jc"i-f3x2-xfi)) ^2 -(101 -2092x^1 -2623X-268 (2^2 -F3;cfi))/9Ho 
-(162 + 138^ + 484/3x2 + 220/3 (2;c"i-xfi))Ho,o- (268 -356;c 
- 88 (5x"^ -2^2 +xf^))/3H2 + (160 -46x + 340/3^2 -778/9 (x"^+xf^)) Hi 
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- 88/3 (l+Jc+Jc2-jcfi)Hi,o + (452- 308x^1 -232x+ 176 (2jc2-xfi))/3Hi,i 
+ 8/3(6+llx-i+6x+llx2)H_i^o-8(l+x)(8^3 + 8C2Ho + 22Ho,o,o-12H3 
+ 4H2,i)+4pgg(-x)(7C3-2C2Ho + 9Ho,o,o-2H3-2H_2,o-8C2H_i+4H_i,2 
-8H_i,_i,o-6H_i,o,o)+4pggW (-31^3 + 2^2 Ho + 6H_2,o-31Ho,o,o + 2H3 
-2C2Hi-20H2,o + 10H2,i+6Hi,o,o + 10Hi,2 + 10Hi,i,o-12Hi,i,i) 
+ 5(1 -;c)(30425/162- 242/3^3 + 830/9^2 + 101/5^1)) 

+ Qn^ ((358x^1 - 928 - 502x + 672xf i)/27 - 1 16/3jc2 + (362 - 260jc~^ - 120jcf ^ 
-238jc+172x2)/9Ho + (36-16x"^-12x + 40(jc2-jcfi))/3Ho,o-(350-168jcf^ 
-88(2x-i+jc-3x2))/9Hi + 8/3((jc-1-1+2jc)(2Hi,o-Hi,i) + (1+x)(H2-^2)) 
-8/3pgg(x)(C2 + 2Hi.o-4Hij+2H2) -5(1 -x)(4112/81 + 8^2 + 28/3^3)) 

+ Cfn^(2;cfi-(744 + 723;c-2017;c2 + 928x"')/27 + (8 + 16/3(x-i+x2) + 12;c)^2 
-8/9(63 + 19x"^ + 81x + 41x2)Ho + (14 + 22jc + 32/3(jc"^+x2))Ho,o-(12 + 8x 
+ 32/3x^1) H2 + (78 -16jc"1-150x + 88x2)/9Hi +4/3 (4x^1 -3x-4jc2 + 3) Hi,i 
+ 4(1+x)(4C3 + 4C2Ho+11Ho,o,o-6H3 + 2H2,i)-5(1-x)(63/2-24C3)) 

+ n^(8/27(l+;c + Jc2-xfi)(5 + 3Hi) + 5(1 -jc)(100/81- 8/9^2)) , (A.3) 

c^2^(jc) = Qn^ (- (2590 + 296x^1 -74x)/27- 695/18x2 + (22 -60x + 56x2) ^2 
+ (16 + 88x)C3 + ((196 + 48x"^-2810x + 582x2)/9+16C2)Ho + (8 + 64x"^ 

+ 212x+116x2)/3Ho,o + (10-32x"i-212x+180x2)/3H2 + (44 + 248x)HoAO 
-8(l+4x + 2x2)H_2,o-(44 + 72x + 40x2)H3-8/3(3 + 2x"^+9x + 8x2)H_i,o 

- (10/3 + 80x^-4/3 (4x^1 +53x)) Hij - (127 + 24x"i - 818x+ 1016x2)/9Hi 
-8Hi,o + 8(4-2x + 6x2)H2,i+4pqg(-x)(^2(3Ho-2H_i)-5H_i,o,o + 2H_i,2) 
+ Pqg(x)(8^2Hi+2/3Hi,o-12H2,o + 16Hi,o,o + 8Hi,2 + 8Hi,i,o-20Hi,i,i)) 

+ CpHf ( + 62 - 71x+ 155/2x2 - (22 - 24x + 4x2) ^2 - (40 - 1 12x + 96x2) ^3 
-(10-47x + 98x2 + 16^2^^)Ho-(7 + 80x-104x2)Ho,o + 4(6-4x+x2)H2 
+ 16 (1 +x) H_i,o + 16 ( 1 + 2x2) H_2,o + (1 + 36x - 48x2) Hi + (6 - 36x + 60x2) ^^^^ 
-4xHij-x2(44Ho,o,o + 8H2,i-24H3)-8pqg(-x)(C2H_i-H_i,o,o + 2H_i,_i,o) 
-Pqg(x)(4C2(Ho + 3Hi)+20Hi,i+22Ho,o,o + 28H2,o-8H3 + 4Hi,o,o-12Hi,2 
-12Hi,i,o + 4Hi,i,i)) 

+ «/((112-392x + 833x2)/27-4/9(34-56x + 53x2)Ho + ;?qg(x)(58/9Hi-4C2 
+ 4/3(Ho,o-H2+Hi,o + Hij))) (A.4) 

As already discussed in Ref. [[T3l . the second moment of cl' +^1' directly enters the NNLO 



Higgs decay rate to hadrons in the heavy-top limit, and agrees with the result of Ref. 11341 . see 
also Ref. [351 . We expect that these coefficient functions will be useful also for other theoretical 
studies. They can be employed, for instance, to extend the large-x results of Ref. Il28l to Pqg . 
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Appendix B 

Since the exact expressions of the NNLO splitting functions are rather lengthy and complex, it 
is useful to have at one's disposal compact but accurate approximate representations which also 
can be transformed readily to Mellin space at all (complex) values of N. In this final appendix 
we therefore provide parametrizations of all NNLO timelike splitting functions in QCD which are 
built up, besides powers of jc, only from the +-distribution and the end-point logarithms 

©0 = 1/(1 -■^)+5 Li=ln(l-x), Lo = Injc. 

('')\'T'-\- 

The non- singlet splitting functions Pns [fT2l| can be represented by 

pi^)^+(;c) = +1174.898©o + 1295.6255(1-x)-707.67Li + 1658.7-4249.4jc 
- 1075.3x2 + 593.9x3 -LoLi[56.907 + 519.37Lo + 559.1x] 
+ 1327. 5Lo- 189.37 L^- 352/9 L^+ 128/81 L^ 
+ /7^(- 183.187 ©0- 173.9355(1 -x)+5120/81Li- 198.1 +466.29X 
+ 181. 18^2 -31. 84^3- 39.1 13xLo-LoLi [50.758 -85.72X 
- 28.55 1 Lo + 23. 102xLo] - 168.89Lo - 176/81 L^ + 64/27 L^ ) 
+ n/ ( -©0- (51/16 + 3^3 -5^2) 5(1 -x)+x(1-x)-1Lo(3/2Lo + 5) 

+ l + (l-jc)(6+ll/2Lo + 3/4Lg))64/81 (B.l) 

and 

Pi^^'^'ix) ^ +1174.898^0 + 1295.6225(1 -x)-707.94Li + 1981.3 -4885.7X 
- 577.42x^ + 407.89x3 +L0L1 [4563.2 + 1905.4Lo - 5140.6x 
+ 1969.5xLo] - 34.683xL^ - 437.03 xL^ + 1625.5Lo - 38.298L^ 
-1024/27^- 140/8 IL^ 
+ n^(- 183.187 ©0- 173.93765(1 -x)+5120/81Li -217. 84 + 511. 92x 
+ 209.19x2- 85.786x3 +92.453xLo + LoLi [71.428 + 30.554 Lo 
-23.722x-18.975xLo]- 188.99 Lo- 784/81 L^+ 128/81 L^) 
+ n/(-©o -(51/16 + 3^3-5^2) 5(1 -x)+x(1-x)-1Lo(3/2Lo + 5) 

+ l + (l-x)(6 + ll/2Lo + 3/4L^))64/81 . (B.2) 

The n ? parts of Pns (which are identical and equal to their spacelike counterparts in Ref. ifTOl '). 
the +-distribution contributions (up to a numerical truncation of the coefficients involving C,i), and 
the rational coefficients of the (sub-)leading regular end-point terms are exact in Eqs. (IB. II) and 
(IB.2I) . The remaining coefficients have been determined by fits to the exact results at 10^^ < x < 
1 — 10^^, and finally the coefficients of 5(1— x) have been adjusted very slightly using the lowest 
integer moments. The difference between the NNLO 'valence' and 'minus' splitting functions is 
equal to that in the spacelike case; a parametrization can be found in Eq. (4.24) of Ref. ifTOl ). 
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Corresponding representations of the pure- singlet and gluon-gluon splitting functions |fT3l are 

P^^^^{x) ^ {n^(-5.926Lf-9.751L?-8.65Li-106.65-848.97x + 368.79x2 
- 61. 2S4 x^ + 96. n 1 LqLi +656 .49 Lo + 425. UlI + 41. 322 lI 
+ 9.012L^ + 419.Slx-^ + 324.01 x-^Lo-12S/9x-^lI-256/9x-^lI) 
+ nj (1.778L? + 16.61 1 U + SI. 195 - 57.688 jc - 41. 827 x^ + 25.628^3 
-7.9934/- 2.1031 L0L1+ 26.294 xLo -7.8645x4 + 57.713Lo 
+ 9.1682L^- 1.9L^ + 0.019122L^- 128/81;c"^)}(l -x) (B.3) 

and 

P^f^ix) ^ +2643.521 ©0 + 4425.451 5(1 -Jc)-3590.1Li -28489 +7469X 

+ 30421 x^- 53017^3 + 19556/ -LoLi(186.4 + 21328Lo) + 12258Lo 
+ 13528 L^ + 3281.7L^ + 191.99 L^ + 5685.8jcL^ + 14214.4x-^ 
+ 10233x"^Lo + 3651.1x"1l^ + 3168x"^L^ + 576x"1l^ 
+ rifi -412.172©o-528.7195(l-x) + 319.97Li+248.95 + 260.6x 
+ 272.79/ + 2133.2/- 926.87/ + L0L1 (1266.5- 29.709 Lo 
+ 87.771Li)+4.9934Lo + 482.94Lg+ 155. 1L^+ 18.085 L^ 
+ 485.18jc4-804.13jc"^-5.47x"% + 2368/9jc"^Lg + 448/9x"^L^) 
+ n/( -16/9©o + 6.46285(l-/-77. 19 + 153. 27X-106.03/+11.995/ 
-LqLi (115.01 -96.522x + 62.908Lo)-69.712Lo-44.8L^-5.037L^ 
+ 472/243 jc"^ + 368/8 lx"^Lo + 32/27 x^^L^) . (B.4) 

In Eqs. (IB.3I) and (IB.4I) the small-x leading terms are exact up to truncations of irrational numbers. 
The same holds for the coefficients of Li, Lj and Lj in Eq. (IB. 31) and that of ©0 in Eq. (IB.4I) . where 
the coefficient of 5( 1 —x) has been minimally adjusted using the lowest moments. 

The new NNLO off-diagonal quantities (fTSl) and (fT9l) can finally be parametrized as 

pjg)^(;c) ^ n^(100/27Lt + 350/9L^ + 263.07L? + 693.84Li + 603.71-882.48x 
+ 4723.2/-4745.8/-175.28/-LoLi (1809.4+ 107. 59jc) 
-885.5xL^+1864Lo + 1512L^ + 361.28L^ + 42.328L^+ 1141.7 jc^i 
+ 615.S3x-^Lo-64x-\lI + lI)) 

+ nj{- 100/27 L^ - 35.446 Lf - 103.609Li - 1 13.81 + 341.26x- 853.35/ 
+ 492.1 / + 14.803/ + L0L1 (966.96- 1.593 Li- 709.1/- 333. 8xL^ 
+ 619.75Lo + 255. 62L^ + 21. 569 L^- 2.8986x^1- 3. 1752x"1Lo 
-32/21 x-^lI) (B.5) 

+ n/(4 + 6(Lo+Li) + (l-2x + 2x2)(3.8696 + 4(Lo + ^i)+3(Lo+Li)2))| 

where the n^ part is exact, and 
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P^l^'^ix) ^ +400/81Lt + 520/27Lf-220.13L?-152.6Li+272.85-7188.7;c 

+ 5693.2x2 +146.98x3 + 128.19/ -LqLi (1300.6 + 71.23Li) + 543. 8;cL^ 
+ 4.4136Lo-0.71252Lg- 126.38L^- 30.061 L^ + 5803.7x-^ 
+ 4776.5x"^Lo + 1001.89x"^L^ + 3712/3x-1l^ + 256x-1l^ 
+ n^(80/81L^ + 1040/81Lf-16.914Li-871.3 + 790.13jc-241.23jc2 
+ 43. 252^3- 4.3465 xL^ + 55.048 LoLi -492 Lo- 343. 1L^-48.6L^ 
+ 6.0041x^1 + 141.93x"% + 2912/27x"1l^ + 1280/81x"1l^) . (B.6) 

The coefficients of the leading small-x terms are exact in Eqs. (IB. 51 ) and (IB. 61 ). up to a truncation 
of irrational numbers. The same holds for the coefficients of Lq, Lq and Lq in Eq. (IB. 61 ). The 
coefficients of the large-x terms are also partially exact. 

Except for values of x very close to zeros of the splitting functions, the parametrizations. (IB. II) - 
(IB. 61) deviate from the exact expressions by less than one part in a thousand, which should be amply 
sufficient for foreseeable numerical applications. Also the complex-A^^ moments of the splitting 
functions can be readily obtained to a perfectly sufficient accuracy using the above representations. 
The Mellin transform of Eqs. (IB. II) - (IB. 61) involve only simple harmonic sums (see, e.g, the 
appendix of Ref. IISSlD of which the analytic continuations in terms of logarithmic derivatives of 
Euler's F-function are well known. 
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